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In this paper we make predictions for nondiagonal parton distributions in 
a proton in the LLA. We calculate the DGLAP-type evolution kernels in the 



| LLA, solve the nondiagonal GLAP evolution equations with a modified version 

of the CTEQ-package and comment on the range of applicability of the LLA 

x ■ 

in the asymmetric regime. We show that the nondiagonal gluon distribution 
g(x\,X2,t, fj?) can be well approximated at small x by the conventional gluon 
density xG(x, fi 2 ). 
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I. INTRODUCTION 

Due to the experimental possibility of probing nondiagonal distributions in hard diffrac- 
tive electro-production processes, theoretical interest in this area in recent years [0-0] has 
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produced interesting results. A pioneering analysis of the nondiagonal distributions for the 
diffractive photoproduction of Z°-bosons in DIS where the applicability of PQCD is guar- 
anteed was given by Bartels and Loewe in 1982 || but went essentially unnoticed. 

In the this paper we would like to complement these results by concrete predictions, albeit 
to the LLA, which can be tested by an experiment. In Sec. |T| we shall demonstrate that in 
the limit of small x the amplitudes of hard diffractive processes can be calculated in terms 
of discontinuities of nondiagonal parton distributions. The real part of the amplitude will 
be calculated by applying a dispersion representation of the amplitude over x. We will show 
that the term in the amplitude which cannot be calculated in terms of the discontinuities 
of nondiagonal parton distributions is suppressed by one power of x in this limit. This 
result will make it possible to calculate the evolution kernels in the LLA following the 
traditional methods |9| and to compare them to results obtained in the QCD-string operator 
approach [ |T0|j . 

In Sec. we calculate the nondiagonal kernels and find them equivalent to those in 
fHH. They are different from the evolution equations for nondiagonal parton densities which 
were presented without derivation in [11]. In Sec. 0we shall make predictions about the 
nondiagonal parton distributions by solving, numerically, the nondiagonal GLAP evolution 
equations with the help of a modified version of the CTEQ-package. In Sec. we shall discuss 
the limitations of the approximation and the need for NLO-results. Future directions will 
be discussed in the conclusions. 



II. NONDIAGONAL PARTON DISTRIBUTIONS AND HARD DIFFRACTIVE 

PROCESSES. 

It has been recently understood that the major difference in QCD between leading twist 
effects in DIS and higher-twist effects in hard diffractive processes is to be attributed to 
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the fact that the latter, initiated by highly virtual, longitudinally polarized photons, can be 
calculated in terms of nondiagonal, rather than diagonal, parton distributions J7|. 

Thus, in order to calculate unambiguously hard two-body processes, it is necessary to 
calculate nondiagonal parton distributions in a nucleon. This implies knowledge of the non- 
perturbative nondiagonal parton distributions in the nucleon which have not been measured 
so far. Hence, the aim of this section is to express the nondiagonal parton distributions in 
the nucleon through quantities being maximally close to the diagonal parton distributions. 
Our second aim is to elucidate on the kinematics of the nondiagonal parton distributions in 
the nucleon needed to describe hard diffractive processes. We shall also discuss the expected, 
limiting, behaviour of the nondiagonal parton distributions. 

For the leading twist effects QCD evolution equations have traditionally been discussed 
in terms of the imaginary part of the amplitude. This is because the bulk of experimental 
data available is on the total cross sections of inclusive processes. However it is well known 
that the QCD evolution equation has a simple form for the whole amplitude which includes 



both real and imaginary parts [12| . This form of the evolution equation can be generalized to 
the case of higher-twist processes, hard diffractive processes || and hard two-body processes 
0]. The analysis of the QCD evolution equation for the nondiagonal parton densities shows 
that the evolution equation contains two terms. The first one is described by a DGLAP- 
type evolution equation [fj-|4],[7]] , whereas the second term, found in Ref. [|J for vector meson 
production at small x, cannot be interpreted in terms of parton distributions. The QCD 
evolution of this term is governed by the Brodsky-Lepage evolution equation 0,0. 

A. GLAP evolution equation for hard diffractive processes 

The aim of this section is to prove that for hard diffractive processes in general, the 
Q 2 -evolution at any x in the DGLAP-region as discussed below, is described by a nondiag- 
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onal GLAP-type evolution equation with asymmetric DGLAP-type kernels and that these 
processes can be calculated through the discontinuity of hard amplitudes. This property is 
important for the quantitative calculations since the dispersive contribution has a relatively 
simple physical interpretation and a deep relation with the conventional parton densities. 
As to the first step, we shall deduce a relationship between amplitudes of hard two-body 
processes and parton densities, and we will find an additional term which has no probabilistic 
interpretation. We will restrict ourselves to the Q 2 -region where the parton distributions are 
still rising and the additional term is of no importance as discussed below. 

The QCD factorization theorem for hard processes means that the hard blob can be 
factorized from the soft one with a precision of a power of Q 2 . The topologically dominant 
Feynman diagrams for small x processes correspond to attachments of only two gluons to the 
hard blob. Although our analysis is rather general, for certainty we shall restrict ourselves 
to the case of diffractive processes where diagrams with two-gluon exchange dominate.^ It is 
convenient to decompose the momentum of the exchanged gluon k in Sudakov-type variables: 

k = xip + (3q + k t , (1) 

where 

p 2 = q 2 = and {k t q) = (k t p) = 0. (2) 

To express the amplitude in terms of non-diagonal parton distributions, the contour of 
integration over j3 should be closed over the singularities of the amplitude in gluon-nucleon 
scattering at fixed x\ and x\ — x. The singularities over ft are located in the complex plane 
of discontinuities over the gluon virtualities: — and and from the s- and u- channel 

° X\ XI— X 1 

discontinuities: tjz^t an d hj^z^ ■ The amplitude differs from if these singularities pinch 
the contour of integration. This causality condition restricts the region of integration to: 



* Hard collisions due to the exchange of 2 quarks are numerically small in the LLA at small x. 
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-l + z<zi<l. (3) 

Our main interest is in the amplitude in the physical region where — t > but small as 
compared to other relevant scales of the process under consideration. In this region the 
amplitude can be represented as the sum of terms having s- or ^-channel singularities only. 
For the s-channel contribution to the amplitude of hard diffractive processes, given 1 > 
X\ > x, the integral over (3 can only be closed over the discontinuity in the amplitude of the 
gluon-nucleon scattering in the variable s. 

Thus this contribution to the amplitude is expressed through the imaginary part of the 
amplitude for gluon-nucleon scattering. The QCD evolution of this term is described by a 
GLAP-type evolution equation where the kernel accounts for the off-diagonal kinematics. 
One also has to add a similar term corresponding to u-channel singularities. 

The contribution of the region x > x± > has no direct relationship to the conventional 
parton densities. This is because the integral over (3 cannot be closed for s- or w-channel 
discontinuities but it may be closed for the discontinuities over the gluon "mass". In Ref. 
the analogy of this term with the wave function of a vector meson has been suggested. 
The presence of this piece which cannot be evaluated in terms of parton densities introduces 
theoretical uncertainties into the treatment of hard two-body processes at large and moderate 
x. 

However, for the imaginary part of the amplitude, more severe restrictions on the region 
of integration apply: 

1 > xi > x. (4) 

In order to ensure that we have not included any superfluous regions of integration, we 
consider the discontinuity of the hard amplitude in the s- channel. This additional restriction 
follows from the requirement that the mass of a qq system in the intermediate state should 
be larger than the mass of a produced vector meson such that the hadron cannot decay into 
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the qq system. This result helps to prove that the piece which cannot be evaluated in terms 
of parton densities is inessential for hard diffractive processes. Let us now apply a dispersion 
representation over the variable s. The only term which cannot be reconstructed in terms 
of a dispersion relation, i.e in terms of discontinuities of parton densities, is the subtraction 
constant^ in the real part. The contribution of the subtraction term to the amplitude with 
a positive signature, i.e symmetric under the transposition of s — > u, is suppressed by 
an additional power of s or, equivalently, by an additional power of x. For the processes 
with negative charge parity in the crossed channel, i.e. electroproduction of a neutral pion, 
the amplitude is antisymmetric under the transposition s — ► wQ. This amplitude has no 
subtraction terms at all, since, in QCD, it increases with energy slower than s^. Thus, in 
this case, a dispersion representation gives the full description. To summarize let us point 
out once more that the small x behaviour of hard diffractive processes is described through 
the discontinuities of hard amplitudes. 



B. Small X{ behaviour of the nondiagonal gluon distribution 

We want to stress that the slope of the x dependence of the amplitudes for diffractive 
processes, however not their residue, should be independent of the asymmetry between frac- 
tions X\ and 22 of the nucleon momentum carried by the initial and final gluons. This is 
due to the fact that the X{ of the partons in the ladder are essential , but not the external 
x, and increase with the length of the parton ladder. Hence, the asymmetry between the 
gluons may be important in one or two rungs of the ladder but not in the whole ladder. 

^This constant is independent of s. 
^This corresponds to a negative signature. 

§An odderon-type contribution in PQCD is suppressed by an additional power of Q 2 . 
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Therefore, at sufficiently small x, it is legitimate to neglect x in most of the rungs of the 
ladder as compared to the X{. This means that the asymmetry between the gluons influences 
the residue but not the slope of the x dependence. 

Let us now discuss the small Xi behaviour of g(x\, x 2 ) - the nondiagonal gluon density 
in a nucleon. The factorization theorem - Eq. 3 of Q- is the basis for the formal definition 
of the nondiagonal gluon density as the matrix element of gauge-invariant bilocal operators 
(cf. eq. 6 of ||): 

g 9/p (x 1 ,x 2 ,t,^) = - r e-W-ir^i TG u + (0,y-,0 T )VG» + (0) \p). (5) 

J -DC 4:71 

Here V is a path ordered exponential of a gluon field along the light-like line joining the 
two gluon operators, t is the square of invariant momentum transferred to the target, and 
\x describes the scale dependence. The sum over transverse gluon polarizations is implied. 
Actually Eq. [| coincides with the definition given in |3|fJ[] for the same quantity. 

For x\ = x 2 , g g / p (xi,x 2 ) is related to the diagonal gluon distribution as xG^^x) — 
^nondiagt^, x). Within the leading a s \iax approximation where the difference between \nxi 
and lnx can be neglected, this distribution coincides with the diagonal one Jj]]: 

^leading a B In x 

(x 1 ,x 2 ,t,fi ) = xG(x,fi ). (6) 

We want to stress here that at fixed Bjorken variable x, the cross sections of hard diffrac- 
tive processes are expressed through g(xi, x 2 , t, /i 2 ) where This can be proved 

by calculating the high energy limit of hard diffractive processes and then applying Ward 
identities similar to Ref. This means that the region of integration near x 2 = (x 2 <ti x±) 
gives only a small contribution to the amplitudes of hard diffractive processes. 

Alternatively, one can examine the leading regions of the integrals in the calculation 
of the distribution of a parton in a parton which is imperative in finding the correct hard 
scattering coefficients for the desired process. This calculation is necessary since one has not 
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only ultraviolet divergences in the partonic cross sections from which one wants to extract the 
Wilson coefficients, but also infrared divergences stemming from initial-state collinearities 



of the participating partons (see Ref. |14j] for further details) which are cancelled by the 
perturbatively calculated expansion of the parton distribution. The region of x 2 = does 
not give a leading contribution. This can be seen by using a simple argument that proper 
Feynman diagrams have no singularity at x 2 = 0, and the region of integration over the 
exchanged gluon momenta x 2 = forms an insignificant part of the permitted phase volume. 

In the first step one has to show that a gluon with x 2 = corresponds to a soft gluon 
and then one can use the argument by Collins and Sterman Jl5| first introduced for proving 
factorization in inclusive e + e~-reactions: 

• For clarification, the quark-loop to which the gluons attach consists both of the hard 
part and the part whose momenta are parallel to the vector meson and of the order 

v _ 2gpjy 

• The minus component Z_ of the quark-loop which is transferred to the target proton 
is /_ « itlnX. The minus component of the gluon momentum k is: 

2qp N 

where m q q is the mass of the qq system. Thus we can neglect k_ with respect to /_ 
in the calculation of the leading term of the amplitude corresponding to the leading 
power in the energy of the process. 

• The transverse momentum in the quark loop (l t ) is cut off by the vector meson wave 
function and thus I 2 -C Q 2 in stark contrast to DIS, whereas the k 2 of the gluon is only 
restricted by the virtualities in the photon wave function which can be as high as Q 2 



8 



0. However one has to satisfy the Ward-identity k^T^ u = where T stands for the 
amplitude. Using Sudakov-variables this becomes: 

x 2 p,T^ + k t ,T^ = 0, (8) 

thus for X2 — > the transverse momentum of the gluon is very small and can be safely 
neglected as compared to l t . 

• One concludes from the above said that the X2 = region corresponds to a soft gluon 
(k 2 ~ 0) and we can use the argument by Collins and Sterman. 

Keeping the above said in mind, the /c-integral involved in the determination of the 
leading regions is of the form (up to overall factors): 

lsoftk dAk {{i-ky + ie){k 2 + i6) f{l ~ 

softk dk+ ((1+ - k + )l. - 1} + ie){2k + k_ - kf + ie) /(Z " k}P)} (9) 



where f(l — k,p) is the amplitude of gluon-nucleon scattering and the integrals over k_ and 
k 2 are suppressed for convenience. If one now integrates over the remaining k + momentum 
one will have the following situations: 

• k + k_ > k 2 : There are no obstructions in the deformation of the integration contour 
since l t k t and therefore the pole at k 2 = is far from the region of integration 
allowed in the LLA. In other words this region does not give a leading contribution. 

** The similarity to DIS will be restored at extremely large Q 2 as a consequence of both a Sudakov- 
type form factor in the photon vertex and a slow decrease, with increasing kt, of the vector meson 
wave function. 
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• k + k_ <C kf: There are no obstructions to the contour deformation since l t 3> k t or in 
other words can be neglected as compared to k\ since the singularities do not 

pinch the contour of integration. Thus this region does not give a leading contribution 
either. 

Thus one has proved that if one of the gluons attaching the soft to the hard part has 
X2 = 0, it will be soft and thus, according to the above reasoning, the = region of 
integration does not give a leading contribution to the parton distribution. 

III. KERNELS IN THE LLA 

There are several possible ways of calculating the evolution kernels to leading order in 
QCD. We first used the traditional approach of calculating the evolution kernels in the 
LLA, via the method of decay cells of e.g., a quark decaying into a quark [0[J2] and using 
cut-diagram techniques to calculate the appropriate Feynman graphs. 

As a cross-check we calculated the first order corrections to the bi-local quark and gluon 
operators on the light-cone, which not only yielded the nondiagonal kernels for the DGLAP 
equation but also the nondiagonal Brodsky-Lepage kernels since we were calculating the 
whole amplitude, not only its imaginary part. However, since we are not interested in those 
kernels at the moment we will not comment on this fact further, let it be said though that 
our results on the Brodsky-Lepage kernels agree with those of Ref . . 

We performed the calculation of the cross-check in a planar gauge i.e q' ■ A = with 
q' 2 ^ 00 and used once more Sudakov-variables: 

" Changes appropriate to the nondiagonal case were made. 

W The advantage of such a physical gauge being that no gluons couple to the operators to first 
order, simplifying the calculations considerably. 
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k = (3p' + aq' + k t , 



(10) 



where 

p' 2 = q' 2 = and (k t q) = (k t p) = 0. (11) 

Since one is neglecting the proton mass one can set p' = P where P is the proton momentum. 

The insertion of the appropriate bi-local operators for quarks and gluons on the light 
cone into the Feynman graphs for first order corrections to those operators, short circuits 
the +-momentum in the graph, which means that the loop variable k has not +-momentum 
(3P but rather X\P (see |14| for more details on calculating one loop corrections to parton 
dsitributions.). This fact eliviates us from the duty of taking the integral over (3. In the 
calculation of the kernels, it remains to take the integral over a which can be done by 
taking the residues and then isolating the leading term multiplying dk\ / k\ and the tree level 
amplitude. This will then yield the kernels in the leading logarithmic approximation.^ 

In the integral over a, one finds three different residues. Two residues stemming from the 
vertical quark or gluon propagators which yield a = ^ and a = ^ giving a contribution 
in the Brodsky-Lepage region, i.e the Brodsky-Lepage kernels, and one stemming from the 
horizontal propagator yielding a = -, — — , which contributes to the DGLAP region, i.e the 
DGLAP kernels. This is analogous to the statements made in Sec. |J. 

After having taken proper care of the definitions of our quark and gluon distributions 
in the amplitudes, we find the following expressions for the nondiagonal evolution kernels, 
where A is given by x\ — x% and corresponds to xbj of e.g., vector- meson production. For 
the quark — > quark transition we get: 

Note that the quark to quark and gluon to gluon kernels also need the self-energy diagrams to 
regulate the two possible collinear singularities, of course, after proper renormalization. 
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The other kernels are computed the same way from the appropriate diagrams: 
[x\ + xi(l — xi) 2 — xfA] 



rv 

P qg ( Xl ,A) = ^N F 

71 



P gq (xi,A) 



Oi, 



C f 



(1-A) 2 
[1 + (1 - Xl f - A] 



7T 



P gg ( Xl , A) = 2N C [- 
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(13) 
(14) 



a?i - A 



21 -x x 2(l-xi)(l-A) 

(15) 



A word conccerning our regularization prescription is in order. In convoluting the above 
kernels, after appropriate scaling of x\ and A with yx, with a nondiagonal parton density, 
one has to replace z\ and Z2 in the regularization integrals with z\ — > (yi — Xi)/yi and 
z 2 — > (z/i — A) /(yi — xi) . This leads to the following regularization prescription as employed 
in the modified version in the CTEQ package in the next section and in agreement with Ref. 



j1 dyi f(yi) f 1 dy l y 1 f(y 1 ) - XifM 

a yi l-xi/y 1+ J Xl y\ V\ — X\ 
nl , (xi-A)f(yx) f 1 dy 1 y 1 f(y 1 ) - xtfixt) 



J xi [yi 



xi)(yt - A) + Jxi y x y\ - x x 
+f(xx) \n f[ ~ ' J 



+ f{x) \n{l- Xl ) 



dyi yifjyi) - Afjxi] 
Xl yi yx - A 



(16) 



(17) 



1 - A, 

For A = one obtains, necessarily, the diagonal kernels, however for the distributions 
q = xxQ(xx,Q 2 ) and g = xxG(xx,Q 2 ), since we chose the definitions of our nondiagonal 
distributions to go into q an g rather then Q and g. We have cross-checked these results 
with those of Ref. via the conversion formulas given by Radyushkin in a recent paper 
0. The formulas given in a recent paper by Ji || do not seem to agree with ours but 
this is only due to a different choice of independent variables used by Ji. After appropriate 
transformations, the formulas of agree with our results pfl. It should be noted however 
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that the kernels from Ref. 0,|5]|| are given for Q and g and not for q and g as we do. One just 
has to multiply the kernels given by those authors for the quark evolution equations with 
Xi/yi after appropriate changes for independent variables of course. Conversion formulas 
between the different notations can be found in Ref. ||. 

The evolution equation for the quantities g(xi,A) and q(x±,A) take the following form 
in our notation: 

dg{x u A,Q 2 ) r 1 dyi 



dliaQ 2 Jxi yi 
dq(xi,A, Q 2 ) f 1 dy x 



f \P gg g(yx, A, Ql) + P gq q( yi , A, Q 2 ) 
P qq q(y u A,Ql)+P qg g( yi ,A,Ql) 



dlnQ 2 J xi yi 

We are interested in the calculation of the asymptotic distribution in terms of the sym- 
metric distribution in the limit of small x and large Q 2 . The reason why this is possible 
is that in this limit the main contribution originates from the non-diagonal distributions 
at Xi,x 2 = x 1 — A with x\ ^> X\. In the case Xi,x 2 ^> A deviations from the diagonal 
distribution are small and can be neglected. 

In the following section we will present the results of our numerical study and show that 
for the case of X\ ^> x 2 — in the kinematic region of practical interest the diagonal and 
nondiagonal distribution will coincide for large Q up to about a factor of 2. 

IV. PREDICTIONS FOR NONDIAGONAL PARTON DISTRIBUTIONS 

Utilizing a modified version of the CTEQ-package, we calculate the evolution of the 
nondiagonal parton distributions, starting from a low Qq=1.6 GeV and with rather flat initial 
distributions for the diagonal and nondiagonal case by using the most recent global CTEQ- 



fit CTEQ4 [18]. The difference between the initial-diagonal and the initial-nondiagonal 
distribution is the factor x multiplying the nondiagonal distribution i.e g(x\, x 2 ) = x\G(xi) in 
the normalization point, in accordance with our earlier argument that the possible difference 
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in the distributions at small x and large Q is only given by the Q 2 - ev °hrtion of g(xi, x 2 , t, ii 2 ). 

We have only considered light quarks, since we are interested in a proton as the initial 
state hadron and the s-quarks are only considered to give a small correction. The following 
figures (see Fig. |1|) show the ratio of the nondiagonal distribution g(x\,x 2 ) to the diagonal 
distribution xG{x) from Q=7 GeV to Q =110 GeV and x 2 from ^ to x\ with x± = 1.1 10~ 4 , 

l.i io~ 3 , i.i io~ 2 . 

The nondiagonal and diagonal distributions agree for x 2 — > x\, i.e. for vanishing asym- 
metry, as expected, and within a deviation of a factor between 0.2 and 1.7, they agree for 
x 2 <C x\. The expectation that there is no In 22, which would give a singluarity for x 2 — > 0, 
is also supported by our numerical calculations. In fact if one takes the x 2 = limit, we find 
firstly (see Fig. [|) that the ratio of nondiagonal to diagonal distribution is finite i.e no In x 2 
infinity and secondly that the evolution of the nondiagonal distribution differs significantly 
in size and shape from the diagonal distribution as first anticipated by Radyushkin in ||. 

Note that at large Q 2 and fixed A C 1 g(xi,x 2 ) is determined by the initial parton 
distributions at xi,x 2 3> A where the validity of the diagonal approximation for g(xi,x 2 ) 
does not depend on our argument in Sec. |T|. The numerical calculation finds that the ratio 



of nondiagonal to diagonal distribution is larger than 1 as anticipated by Radyushkin ||T7 
based on general arguments about the nature of the double distribution which he discusses 
in§. 

To see whether our numbers i.e. our numerical methods could be trusted, we used the 
MATHEMATICA program to calculate the first iteration and the first derivative of the 
evolution to see how good or bad our numbers were. As it turns out our integration routines 
produce a very good agreement with the numbers from MATHEMATICA with a relative 
difference of 5%. This leads us to believe that our numbers can be trusted to high accuracy 
for x 2 of 0(x\) and within 5% at x 2 down by two orders of magnitude as compared to X\. 

A few words about the nature of the modifications to the CTEQ-package are in order at 
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this point. The basic idea we employed was the following: In the CTEQ package the parton 
distributions are given on a dynamical x- and Q-grid of variable size where the convolution of 
the kernels with the initial distribution is performed on the x-grid. Due to the possibility of 
singular behaviour of the integrands, we perform the convolution integrals by first splitting 
up the region of integration according to the number of grid points, analytically integrating 
between two grid points Xi and Xi + \ and then adding up the contributions from the small 
intervalls. We can do the integration analytically between two neighbouring gridpoints by 
approximating the distribution function through a second order polynomial ay 2 + by + c, 
using the fact that we know the function on the gridpoints Xi-\,Xi and Xi+i and can thus 
compute the coefficients a,b,c of the polynomial. This approximation is warranted if the 
function is well behaved and the neighbouring gridpoints are close together. We treat the 
last integration between the points X\ and Xi (which are not to be confused with the X\ and 
x 2 of the parton ladder) by taking the average of x\ and Xi and the values of the function 
at x\ and 22 and using those together with xi, X2 and the value of the function at x\ and 
£2 to compute the coefficients of the polynomial. The coefficients are computed in the new 
subroutine NEWARRAY and the integration of the different terms in the kernels is performed 
in the new subroutine NINTEGR. The case x\ = A = x is implemented analytically but 
separately in NINTEGR. Appropriate changes in the subroutines NSRHSM, NSRHSP and 
SNRHS were made to accomodate the fact that the kernels and also the integration routines 
changed from the original CTEQ package. A detailed description of the code will be provided 
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V. LIMITATIONS OF THE LLA IN THE NONDIAGONAL CASE 

The LLA approach of the previous sections accounts for the contribution of a certain 
rather limited range of integration in the parton distributions. Regions outside these limits 
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might contribute to the leading power. Looking at some other physical quantities such as F 2 , 
where one finds substantial modifications due to the NLO-terms, we are forced to assume 
that this may be also true in our case. This results in the urgent need to carry out a NLO 
calculation and numerical study of the evolution equation, which will be the next step of our 
program. 



VI. CONCLUSIONS AND OUTLOOK 

In summary, we have calculated the evolution kernels for non-diagonal parton distribu- 
tions in the LLA using traditional methods and found agreement with the results of [|3||5|,|T0 



deduced by other methods. It was important to show that the traditional approach can still 
be applied. Thus traditional methods can be used to calculate systematically hard diffractive 
processes within the NLO approximations. We have also proved the similarity between the 
diagonal and nondiagonal parton distributions. The latter ones determine the cross sections 
of hard diffractive processes in the small x region. We have made predictions about the 
nondiagonal parton distributions within the LLA with the help of a modified version of the 
CTEQ-package. Numerical calculations found the diagonal and nondiagonal gluon distribu- 
tions, which dominate hard diffractive processes, to be very similar at small x as expected 
from the previous discussion. 
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FIG. 1. The fraction g{x\, X2) / x\G{xi) as a function of X2 for fixed xi and various energies Q. 
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FIG. 2. The fraction g(x±, X2) / x\G{x\) as a function of x\ for X2 = and various energies Q. 
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